Zeros of the Potts Model Partition Function in the Large-g Limit 
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We study the zeros of the q-state Potts model partition function Z{A, q, v) for large q, where v 
is the temperature variable and A is a section of a regular d-dimensional lattice with coordination 
number ka and various boundary conditions. We consider the simultaneous thermodynamic limit 
and q ^ oo limit and show that when these limits are taken appropriately, the zeros lie on the unit 



circle |a;A| = 1 in the complex xa plane, where xa 



vq 



-2/k.a 



For large finite sections of some 



lattices we also determine the circular loci near which the zeros lie for large q. 
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PACS numbers: 05.20.-y, 64.60.Cn, 75.10.Hk 



I. INTRODUCTION 

In this paper we shall study the g-state Potts model in 
the limit q —^ oo and shall present some new results on 
the complex-temperature zeros of the partition function 
in this limit. The Potts model has served as a valuable 
model for the study of phase transitions and critical phe- 
nomena P, 01 • The q ^ oo limit of the model is exactly 
solvable on any lattice 0], as can be seen from the fact 
that in this limit the model essentially reduces to a single- 
site problem. Since the q —^ oo limit of the Potts model 
is solvable, it is natural to investigate the properties of 
the model for large q. Indeed, although we will focus 
on partition function zeros rather than physical applica- 
tions, we note that the Potts model with large q has been 
used in modelling the kinetic behavior of soap froths . 

The (zero-field) Potts model is defined, for tempera- 
ture T on a lattice A, or more generally, a graph G, by 
the partition function 



Z{G,q,v) = exp(A'^(5^,^J 



(1.1) 



where Ci — 1 , . . . , g are the classical spin variables on 
each vertex (site) i £ G, {ij) denotes pairs of adjacent 
vertices, K = (3J where /3 = {kBT)~^, and J is the 
spin-spin coupling. We use the notation a ~ , and 
V = a — 1 ^ e'^ — 1. For the Potts ferromagnet, J > 
so w > 0, while for the antiferromagnet, J < so that 
— 1 < u < 0. The graph G is formally defined by its set 
of vertices V and its set of edges (= bonds) E; we de- 
note the number of vertices oi G a.s n — n{G) — \V\ and 
the number of edges of G as e(G) ~ \E\. The complex- 
temperature zeros of Z{G, q, v), i.e., the zeros in the com- 
plex V plane, have the property that on a regular lattice 
graph with dimensionality greater than the lower critical 
dimensionality, d > 2, in the thermodynamic limit Lj — s- 



oo with Lj/Lk a finite nonzero constant, 1 < j, k < d, 
they merge to form a locus B consisting of curves that 
bound regions of different analytic behavior for the cor- 
responding free energy F = —ksTf, where the dimcn- 



sionless function / is given by / = lim„_ 



InZ. For 



the thermodynamic limit of a graph which is a section 
of a regular lattice graph A with some boundary condi- 
tions, we denote ka as the coordination number satisfying 



KA = 2 lim. 



(G)- 



e(G)/n(G) in this limit (independent 



of boundary conditions). Reviews of the Potts model 
include 0,13. 

Let G' = {V,E') be a spanning subgraph of G, i.e. a 
subgraph having the same vertex set V and a subset of 
the edge set, E' C E. Then a useful representation of 
Z{G,q,v) is [| 



Z{G,q,v) 



G'CG 



qk{G')^e{G') 



(1.2) 



where k(G') denotes the number of connected compo- 
nents of G'. Since we only consider connected lattice 
graphs G, we have k{G) = 1. For the ferromagnet, eq. 
(|1.2f) allows one to generalize q from the positive integers 
to the positive reals, K+. As is evident from eq. (|1.2|l . 
Z{G, q, v) is a polynomial in q with minimal and maximal 
degrees 1 and n{G), and a polynomial in v with minimal 
and maximal degrees and e(G). Thus, for a finite G, 
Z{G, q, v) is completely determined by its zeros and can 

be written as Z{G,q,v) = OilTb ~ ^i(9)]- 

We first recall some previous related work. For the 
special case g = 2 on the square lattice, where one has 
the exact Onsager solution for the partition function, it 
was found that (in the thermodynamic limit) the locus B 
consists of the union of circles |(i± 1| — \/2 @. The locus 
B has also been determined exactly for g = 2 on various 
other d = 2 lattices (e.g., [3,I3)- However, aside from 
the special case g = 2, d = 2, one does not, in general. 



know B for lattice graphs of dimensionality d > 2. In- 
deed, studies of complex-temperature zeros of the Potts 
model partition function for finite sections of various lat- 
tices have found that they exhibit a considerable scatter, 
especially in the half-plane with Re{a) < 0-^3 • ^ 
complementary approach has been to obtain exact solu- 
tions for the partition function and finite-width strips of 
lattices with various boundary conditions and hence de- 
termine exactly the locus B in the infinite-length limit 
(e.g., 0|-|2l|- These systems are quasi-one-dimensional, 
so that B does not cross the positive real a axis, i.e., the 
free energy is analytic for all nonzero temperatures. With 
the definition Xsq = v/^/q for the square lattice, it was 
found the complex-temperature zeros calculated for fi- 
nite sections of this square lattice with duality-preserving 
boundary conditions have the property that a subset lies 
exactly on an arc of the unit circle \xsq \ = 1 in the com- 
plex Xsq plane for Re{xsq) ^0^^^,'22]. For phys- 
ical temperature the coefficients of powers of a = are 
positive, so there are no zeros on the positive real axis 
Re{a) > for any finite lattice. However, in the thermo- 
dynamic limit, the phase boundary on B crosses this axis 
at the ferromagnetic phase transition point, ac — 1 + ^fq^ 
i.e. Xc = 1. A number of interesting results have been 
obtained by F. Y. Wu and collaborators: in Ref. ^llj a 
conjecture was made that for finite sections of the square 
lattice with self-dual boundary conditions and for the 
same lattice with free or periodic boundary conditions 
in the thermodynamic limit, the complex-temperature 
zeros of the partition function in the Re{xsq) > half- 
plane are located on the unit circle \xsq\ — 1- In Ref. 
[I2|, it was shown, using Eulcr's identity for integer par- 
titions, that for the Potts model on the self-dual square 
lattice, if one takes the limit g ^ oo, all of the zeros 
of the partition function are located on this unit circle 
\xsq\ — 1 for any lattice size. In Refs. ^2, 13 113 it was 
shown that partition function for the Potts model on a 
finite c?-dimensional Cartesian lattice E'' with d > 2 and 
appropriate boundary conditions, in the limit (7 — > cx), 
divided by an appropriate power of q, is given by the gen- 
erating function of restricted integer partitions in d — 1 
dimensions, with expansion variable t ~ v'^ /q. 

In the present paper we present some results on the 
zeros of the Potts model partition function for large q 
on general regular d-dimcnsional lattices, both in the 
thermodynamic limit and for finite lattices, with various 
boundary conditions. We shall give a simple proof for an 
arbitrary regular lattice of dimension d > 2 that, when 
the thermodynamic limit and the q — > 00 limit are taken 
in an appropriate way, the zeros of the partition function 
lie on the circle \xa\ = 1. This generalizes the corre- 
sponding result of Wu and coworkers from a Cartesian 
lattice E'' to an arbitrary regular d-dimensional lattice. 
We also consider finite sections of some regular lattices 
and determine the circular loci near which the zeros lie 
for large q. We shall often call this the "approximating 
circle" for the zeros. 



II. GENERAL RESULTS FOR LARGE q 

In this section we analyze the pattern of partition 
function zeros of the Potts model for large q. Let us 
first arrange the terms in eq. (|1.2fl in order of de- 
creasing power of V, from ti'^('^) to 1. There are ('^^^■*) 

terms in the coefficients of and , so that in 

the special case 9=1 the partition function reduces to 
Z{G,q = l,v) ^ {v + l)^(^). These {^f>) terms can 
have different powers in q except for very small or very 
large j, as follows. The term with the highest power of 
V arises from the contribution of G" = G in eq. (|1.2|l . 
where all of the edges of G are present, and the power 
of V is e(G); since this graph consists of a single con- 
nected component, the coefficient is q and the term is 
g^e(G)^ In considering other terms in Z{G, q,v), we as- 
sume that the graph G has the property that the removal 
of a few edges does not cut it into disconnected com- 
ponents. This excludes the cases of a one-dimensional 
lattice and quasi-one-dimensional lattice strips of narrow 
widths; we comment on these in the appendix. For exam- 
ple, for a regular c?-dimensional lattice graph, we require 
that d > 2 and that the lengths Lk with 1 < fc < c? are 
great enough so that this condition is satisfied. Next, 
the term in Z{G,q,v) with one lower power of v is ob- 
tained from the set of G' each of which has an edge set E' 
with one less edge than the full edge set E. This term is 
("^'f)) where (^) = ^!/[m!(€ - m)!]. Similarly, 

the first few terms with descending powers of v below 
the maximal power are given by ('^'^■') qv'^^^^^^ , where 
j edges are removed in G' and the linear dependence on 
q reflects the property that this edge removal maintains 
the connectedness of the resulting G', so fc(G') = 1. In 
general, this holds for j < <5(G), where 6{G) and A(G) 
denote the minimal and maximal degrees of a vertex in 
G. (The degree of a vertex is defined as the number of 
edges connected to it). Thus, for example, a term in- 
volving q^ begins to appear in the coefficient of 
for cyclic strips of the square lattice with L\ = > 2 
and L2 = Ly > 2 since the degree of vertices on the 
upper and lower boundaries is three, but it appears in 
the coefficient of for the same strip with = 2. 

At the other end of the polynomial Z{G, q, v) ordered 
as powers of v, the term with the power of v equal to 
zero arises from the contribution of G' with E' — 0, i.e., 
the spanning subgraph containing no edges, leaving the 
n(G) vertices as disconnected components; this term is 
thus q^^'^'^K The term linear in v results from the contri- 
butions of spanning subgraphs G' with one edge and is 
e{G) q^^'^^~'^ V. In general, the first few such terms are 
given by {^'^^^^ q"^^^~^ , where G' consist of j edges 
such that j < g{G), with g{G) denoting the girth of G, 
i.e., the length of minimum closed circuit on G. For ex- 
ample, for cyclic strips of the square lattice again, the 
term appears first in the coefficient of for 

Lj; > 3 because the girth is four, but term ap- 
pears first in the coefficient of if La, = 3 and the term 



2 



qn(G}-i appears in the coefficient of if = 2. Thus, 
for graphs G satisfying the above-mentioned conditions, 
we have 

e(G) 

(2.1) 

where the coefficients Cj{q) are polynomials in q. 

Now let us consider the limit q ^ oo. Evidently, in 
this limit, for a given finite graph G, there is a single 
dominant term in Z(G, q, v), namely g", arising from G" 
with E' = %. This term is what one would get in the 
evaluation of Z if there were no spin-spin interactions, 
but instead just single-site contributions. This reduction 
helps to understand the fact that the Potts model is ex- 
actly solvable '3| in the limit q ^ oo. 

Next, let us examine the terms in Z{G, q, v) in more 
detail. For a given graph G, we define 



XG 



MG)/e{G) 



and, in the thermodynamic limit, 

V 

„Z/KA 



(2.2) 



(2.3) 



Considering q to be large, we focus on the first two and 
last two of the terms in eq. (|2.1|l : we can write these as 



Z{G,q,v) 



= qn{G) + l 



4(^){l-Fe(G)[xGg"(^^/^(^)]-i} + 



(2.4) 



We can approximate Z{G, q, v) by including just the first 
and last terms if several conditions are met. In the ex- 
pression multiplying Xq^\ the second term is negligible 
compared with the first if xq ^ 0(1) and the condition 



q » e(G)^('^)/("(«)-i) 



(2.5) 



is satisfied. In the expression multiplying q^^, the first 
term is negligible compared with the second if xg ^ 0(1) 
and the condition 



q > [e(G)]''('^)/(^('^)~"(<^)+i) 



(2.6) 



is satisfied. Let us denote the degree of Cj{q) in q as pj. 
Among the ('^^^^) terms in Cj{q), only a few have this 
power except for small and large j as discussed above. 



Now we require that Cj{q)v^ be small compared with 5", 
so that in terms of the power of q, 



or 



Pj < 



e(G) 

(e(G)-j)(n(G)-l) 
e(G) 



(2.7) 



(2., 



Further, if these conditions hold and xq ^ 0(1), then 
terms arising from the third, fourth, and subsequent 
terms of Z(G,q,v) are small relative to the first, and 
the third, fourth, and earlier terms counting in from the 
last, are negligible relative to the last, with the ordering 
specified by ascending powers of w, as in eq. H2.1|l . Pro- 
vided that Xg ^ 0(1) and the above conditions hold, one 
can approximate Z{G, q, v) by keeping only the first and 
last terms in ec^. (I2.1|) : 

Z(G,(z,t;)-(?w"(^)+g"(^) = q"^^'>+^[x'(^^^ +q-^] . (2.9) 
The zeros of this approximation to Z{G,q,v) are given 

by 



XG = (~9 ) 



lNl/e(G) 



(2.10) 



Now consider the limit e(G) — > 00 and n{G) 00 with 
the ratio e{G)/n{G) = kg/2 finite. If q also goes to 
infinity, sufficiently fast to satisfy the conditions H2.5I) 
and (Eini) but more slowly than the exponential e ' 
with b a real positive constant, then 



lim lim | q 

q-^oo e(G)~>oo 



l/e(G) 



and these zeros merge onto the unit circle 



\xg\ 



1 



(2.11) 



(2.12) 



If q were to grow more rapidly, as g ^ ^^'^i'^) -where h is 
an arbitrary positive real constant, then the zeros would 
merge onto the circle \xg\ = e~^. 

In particular, one may consider the thermodynamic 
limit of a graph G which is a section of the regular lattice 
A. Then the conditions corresponding to eqs. (|2.5(l and 
(|2.t)|) are, respectively. 



ka/2 



and 



q > e(G) 



q > e(G)'''^/('"'^"2) 



(2.13) 



(2.14) 



The inequality ka/2 > ka/(ka — 2) is equivalent to the 
inequality ka > 4. Since ka > 4 for all of the lattices 
that we consider except for the honeycomb lattice, a con- 
sequence is that for these former lattices with k\ > 4, 
condition (|2.13|) implies condition (|2.14|) . Provided that 
x\ 0(1) and the above conditions hold, so that one 
can approximate Z{G,q,v) with the first and last terms 



3 



in eq. (|2.1() . it follows that in this thermodynamic limit, 
with q also going to infinity sufficiently fast to satisfy the 
conditions H2.13|l and (|2.14|l but more slowly than the ex- 
ponential form given above, the zeros of Z{G, q, v) merge 
onto the circle 



\xa\ = 1 



(2.15) 



where xa is given by eq. Our result (|TT5|) holds for 

any regular lattice with dimension d > 2. For example, 
for the d-dimensional Cartestian lattice E'', Kj;d = 2d, 
while for the d-dimensional body-centered cubic lattice. 



lib. 



Since the free energy is nonanalytic at xj^ 
(12.15(1 yields the asymptotic relations 



Vc,A 



as q 



and hence 



■ In g as g — > oo 



1, eq. 



(2.16) 



(2.17) 



for the values of v and K where the Potts model on 
the lattice A has a phase transition from a paramagnetic 
(PM) high-temperature phase to a ferromagnetic (FM) 
low-temperature phase in the limit of large q. We note 
that this agrees with the large-g limit of the mean-field 
theory result [2fil | 



K, 



2(g-l) , . 

c,MFT = — 7 in^g 

KA(q - 2) 



1) 



(2.18) 



VIZ., 



c,MFT 



{l + q-^ +0{q-'^)}\nq 



{q-' + 0{q-^)} 



(2.19) 



as g oo. In cqs. I|2.17|l and l|2.18() . the property that 
K(. increases with increasing q can be understood as a 
consequence of the fact that as q gets large, each spin 
has more possible values (is "floppier"), and hence one 
must go to a lower temperature for the ferromagnetic 
long-range order to occur. The feature that Kc increases 
asymptotically like Ing as g — > oo can be understood 
since an order-disorder transition involves a balance be- 
tween minimizing the configurational energy and max- 
imizing the entropy terms in the free energy per site, 
F = U — TS, and in this limit, the entropy per site is 
S kB^aq. 

Next, we recall the known exact equations for the PM- 
FM phase transition on the 2D lattices 1, 2], which we 
write in a convenient manner for the discussion of the 
large-g limit. 



4-1 

V 



for A = sq 
for A = tri 



1 



3g 



for A — he , 



(2.20) 



where sq, tri, and he denote the square, triangular, and 
honeycomb lattices. As is evident, in the limit of large g 



and V such that g 



for the square lattice, g 



for 



the triangular lattice and g = v^^^ for the honeycomb lat- 
tice, these equations are all in agreement with the general 
relation H2.16|l . which also holds for higher-dimensional 
lattices. 



III. PARTITION FUNCTION ZEROS FOR 
FINITE LATTICE SECTIONS 

A. General Structure 

It is also of interest to consider zeros of Z{G,q,v) on 
finite lattice graphs in the limit of large g. In general, if 
the conditions specified in the previous section are satis- 
fied, one can approximate Z{G, g, v) by keeping the first 
and last terms in eq. (|2.1I) . We have studied these ze- 
ros in detail for sections of various 2D lattices and have 
found that a sufficient criterion that these conditions can 
be satisfied for these and higher-dimensional lattices is 
that one uses fully periodic boundary conditions. Ac- 
cordingly, when g is large, the zeros in the xa plane are 
located close to a circle. For finite lattice graphs with 
fully periodic boundary conditions, e{G)/n{G) — ka/2so 
that the two definitions in eas. (|2.2() and (|2.3I) are equiv- 
alent. However, for finite graphs with other boundary 
conditions, xq varies from one graph to another. There- 
fore, we will plot zeros in the xa plane for the lattice A 
with different boundary conditions. We first rewrite the 
approximate expression (|2.9|l in an equivalent form that 
will be convenient for our discussion; when g is large, 

Z{G,q,v) - gw'^^'^) + g"(^) 

^ ^(2e(G)/KA) + l[^e(G) _^ ^n(G)-(2e(G)/«:A)-l] (g j) 

The zeros of this approximation to Z[G, g, v) are located 
on a circle in the xa plane with radius 



r(A,SC) = gP(^'^^) 



where 



piA,BC) 



n{G) - 1 2 
e{G) KA 



(3.2) 



(3.3) 



(This power p should not be confused with the powers pj 
in eqs. l(T7|) . lfTH|) . The radius r{A,BC) is less (greater) 
than unity if p{A,BC) is negative (positive). For lat- 
tices with periodic boundary conditions in all directions, 

n{G)/e{G) = 2/ka and hence 



piA,PBC) 



1 



e(G) 



(3.4) 
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as in eq. H2.10|l . 

In our explicit calculations of zeros for various lattices, 
we have found that for large but finite q, the approximat- 
ing circles near which these zeros lie may have centers 
that are slightly shifted to the left or right of the origin 
of the xa plane and move in toward the origin as (7 — > oo. 
The following analysis provides an understanding of these 
shifts. As discussed in Sectional there are S{G) terms in 



the partition function given by 



qv 



•'^'^^ ^ and there 



are g{G) terms given by ^v^. For the fi- 

nite lattice graphs with large vertex degree and small 
girth, a specific approximation for the partition function 
Z{G, q, v) when q is large is 



Z{G,q,v) - q{v + ly^^^ + q''^^'> 



(3.5) 



so that the center of the circle, denoted by c{A,BC), is 
not at the origin but at 



c{A,BG) 



1 



(3.6) 



Of course, this approaches zero as g — s- cx), in agreement 
with our result above, eq. (|2.15|l . The triangular lattice 
is an example for which the center of the circle is shifted 
to the left, as in eq. (|3.6|1 and will be discussed below. 

On the other hand, for finite lattice graphs with small 
vertex degree and large girth, a specific approximation 
for Z{G, q, v) when q is large is 

Z(G, q, v) - gw^(^) + q^(G)-e(G) ^ ^y{G) 
= g"(G) + l[^e(G) ^ ^ q{n(G)/e(G))-l^^yiG)^ 

(3.7) 

The circle near which the zeros lie is now shifted to the 
right. Let us first consider the case with fully periodic 
boundary conditions so that xg = xa- In the complex 
xa plane, the circle has radius 



r{K,PBC) = 
and center at 

c(A, PEC) 



,l/e(G) 



q2/e{,G) _ g(2n(G)/e(G))-2 



„n/e(G)-l 



^2/e(G) _ ^(2n(G)/e(G))-2 



(3.8) 



(3.9) 



For the lattices with e(G) larger than n{G), the term 
g(2n(G)/e(G))-2 jj-^ ^]-^g denominator is negligible when q 
is large, and the radius of the circle is approximately 
as in eq. H3.2|l . For large q, the center of the 



q 



-l/e(G) 



circle in the xa plane is 

c(A, PEG) = q-'(^'BC) 

where 

2 - 71(G) 



s{A,PEG) = 1 



e(G) 



(3.10) 
(3.11) 



When e(G) is large, so that 2/e(G) ^ 1, the position 
of center can be further approximated as c(A, PEG) ~ 
^(2/ka)-i^ which approaches the origin as (7 — > 00. 

For the lattice sections with other boundary condi- 
tions, a specific approximation for Z{G,q,v) for large q 
is 

Z(G,g,«)^9^2e(G)/KA) + lx 

[^e(G) ^ ^„(G)-(2e(G)/KA)-l(i + qiV ^^yiG)^ 

(3.12) 

In the Xa plane, the resultant approximating circle has 
radius 



r{A,EG) = 
and center 
c{A,BG) = 



^l/e(G)-(n(G)/e(G))+2/KA 



g2/e(G)-(2n(G)/e(G)) + (4/«:A) _ g(4/KA)-2 



(3.13) 



(2/ka)-1 



^2/e(G)-(2n(G)/e(G)) + (4/KA) _ q(4/KA)-2 ' 

(3.14) 

For lattices with ka > 2, q'*/'^'^ ^ in the denominator is 
negligible when q is large, and the radius of the circle is 
the same as eqs. (13.211 and (|3.3|) . We also have 



s{A,EG) = 1 



2 - 2n(G) 



2 

KA 



— (3.15) 



for large q. The honeycomb lattice is an example for 
which the circle near which the zeros lie shifts to the 
right of the origin and will be discussed below. An- 
other example is the square lattice with next-nearest- 
neighbor interactions that have the same strength as the 
nearest-neighbor interactions, which we have also ana- 
lyzed. From the above discussion, we expect that the 
center of the circle is the origin for lattice sections with 
vertex degree roughly equal to girth. 



B. Square Lattice Sections 

For the square lattice with toroidal (tor) boundary con- 
ditions, e{sq,tor) — 2LxLy, Ksq — 4, and 



p{sq, tor) 



1 



ILxLy 



(3.16) 



Hence the radius of the approximating circle increases as 
the area of the lattice increases, and approaches unity 
from below in the thermodynamic limit. The zeros of 
the partition function for the Potts model on the toroidal 
strip of the square lattice with Ly = A and L^: = 9 in the 
Xsq plane when q = 1000 are shown in Fig. ^a). The 
radius of the circle is about 0.9085. The zeros lie very 
close to this circle for the above value of and we find 
that they get closer to the circle given by eqs. (|3.2|l and 
(|3.3|) when q increases. 
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(a) 




(b) 












(c) 




(d) 



FIG. 1: Zeros of the Potts model partition function, plotted in 
the Xsq plane, for a section of the square lattice with Ly = A and 
Lx = 9, when q = 1000, for the following boundary conditions: (a) 
toroidal; (b) cyclic; (c) cylindrical; (d) free. For comparison, the 
unit circle and circles where the zeros cluster are also shown. 



For the square lattice with cycHc boundary conditio 
we have e{sq, eye) = Lx{2Ly — 1). By ea. H3.3|l . 



p{sq, eye) 



Lx 



2Lx{2Ly - 1) 



(3.17) 



In the thermodynamic hmit the radius of the approxi- 
mating circle decreases and approaches unity from above. 
Note that if one were to keep Ly fixed and take oo, 
thereby violating the premises of our analysis, the radius 
of this approximating circle would be gi/^C^^-y-i)^ which 
would diverge as q oo. The zeros of the cyclic strip 
with Ly = 4 and L^, = 9 in the Xgq plane when q = 1000 
are shown in Fig. Cfb) . The radius of the approximat- 
ing circle is about 1.47. The zeros are not as close to 
the circle as for the corresponding toroidal case. On the 
other hand, if one keeps L^ fixed and takes — > oo, 
corresponding to cylindrical boundary conditions, the ra- 
dius of the approximating circle approaches unity from 
above. The zeros of the cylindrical strip with Ly — A 
and Lj. — 9 in the x^q plane when q — 1000 are shown 
in Fig. n^c). The radius of this approximating circle is 
about 1.11. The zeros are closer to the circle than the 
cyclic strip with the same Ly and Lj. because the num- 
ber of boundary vertices (i.e. vertices with coordination 
number three) is reduced. 

For the square lattice with free boundary conditions, 
e(sg, free) = 2LxLy — Lx — Ly. By ea. (|3.3|l . 



p{sqjree) 



Lx 



Ly -2 



2(2LxLy — Lx 



Ly) 



(3.18) 



In the thermodynamic limit the radius of the circle de- 
creases and approaches unity from above. Again, how- 
ever, if one were to keep Ly fixed and take Lx oo, the 
radius of the approximating circle would be q'^/^('^^v~^\ 
(the same as the cyclic case), which would diverge as 
g — > oo. Since for a given lattice size the number of 
boundary vertices is larger than for other boundary con- 
ditions, and especially since there are four vertices with 
coordination number equal to two, it is natural to expect 
that the zeros will lie farther from the asymptotic circle 
than for these other boundary conditions. The zeros of 
the free strip with Ly = A and L^; = 9 in the Xgq plane 
when q = 1000 are shown in Fig. Q^d). The radius of the 
approximating circle is about 1.90. 



C. Triangular lattice 

For the triangular lattice with toroidal boundary con- 
ditions, e{tri,tor) = SLxLy and K{tri,tor) = 6. Hence, 
for a finite section of the triangular lattice with toroidal 
boundary conditions, as for the infinite triangular lattice, 

V 

:t/3 



Xtr 



q 



(3.19) 



and 



p{tri, tor) 



'dLxLy 



(3.20) 



as given by ea. (|2.10|l . This radius increases as the size 
of the lattice increases, and approaches unity from below 
in the thermodynamic limit. There is a noticeable shift 
of the circle to the left. The center of the approximating 
circle is at Xtri = —q~^^^, approaching zero as q — > oo. 
The zeros of the toroidal strip with Ly — 3 and ^ 9 
in the xtri plane when q — 1000 are shown in Fig. I2a). 
The radius of this approximating circle is about 0.918. 

For a section of the triangular lattice with cyclic 
boundary conditions, e{tri, eye) = L^i^Ly — 2). When q 
is large, we calculate 

for the approximating circle. In the thermodynamic 
limit, the radius of this circle decreases, approaching 
unity from above. However, if one were to keep Ly fixed 
and take ^ co, then the radius would be g2/3(3L„-2)^ 
which would diverge as g ^ oo. The zeros of the cyclic 
strip with Ly = A and La, = 9 in the Xtri plane when 
q = 1000 are shown in Fig. [2Jb). The radius of this ap- 
proximating circle is about 1.47 as for the corresponding 
case of the square lattice. If one keeps fixed and takes 
Ly ^ oo, corresponding to cylindrical boundary condi- 
tions, the radius of the approximating circle approaches 
unity from above. The zeros of the cylindrical strip with 
Lj, = 4 and L^, = 9 in the xtri plane when q = 1000 are 
shown in Fig. ^c). The radius of this approximating 
circle is about 1.12. 

In general, the approximation in ea. l|2.9|l is not valid 
for the triangular lattice with free boundary conditions. 
The reason is that there are two vertices with degree 
two, so that the coefficient of ^e(tri,/ree)-2 contains the 
term 2q^ and the coefficient of ye(trijree)-4 contains 
the term for Lx,Ly > 2. Therefore, when q is 
large there are two pair of roots around v ~ a^/q, or 
equivalently around ±z gi/s in the Xtri plane. We have 
e{tri, free) = SL^Ly — 2Lx — 2Ly + 1. The partition 
function zeros for the free strip of the triangular lattice 
with Ly = 4: and L^ = 9 when q = 1000 are shown in Fig. 
|2td). Although the approximation (|2.9|) does not apply 
here, one can see that many of these zeros lie close to an 
approximating circle with radius 1.84, which is similar 
to what one would get by formally using p{tri, free) ~ 
2{L^ + Ly - 2)/[3{3L^Ly - 2L^ - 2Ly + 1)]. We observe 
that there are two almost overlapping complex-conjugate 
pairs of zeros which have relatively large magnitude and 
are close to the imaginary axis. 

D. Honeycomb lattice 

For the honeycomb lattice with toroidal boundary con- 
ditions where both Lj, and Ly are even, e{he,tor) — 
3LxLy/2 so that n{hc,tor)/e{he,tor) — 2/3 for any L,j; 
and Ly. Thus, 

Xhc = -Z^ (3.22) 




(b) 



(d) 

FIG. 2: Zeros of the Potts model partition function in the Xtri 
plane for the triangular lattice when q = 1000, for the following 
boundary conditions and sizes: (a) toroidal with Ly = 3 and Lx = 
J 9; (b) cyclic with Ly = 4 and Lx = 9; (c) cylindrical with Ly = A 
and Lx = 9; (d) free with Ly = A and Lx = 9. For comparison, 
the unit circle and approximating circles arc also shown. 
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(3.23) 



as given in eg. (|2. 10(1 . This radius increases as the size of 
the lattice increases, and approaches unity from below in 
the thermodynamic limit. There is a shift of the circle to 
the right. By ea. (|3.11|) . 
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1 




(3.24) 



which approaches 1/3 as the size of the lattice section 
becomes large. The zeros of the toroidal strip with Ly = 
4 and m = 9{Lx = 18) in the Xhc plane when q — 1000 
are shown in Fig. Ela). The radius of this approximating 
circle is about 0.938. 

For a section of the honeycomb lattice with cyclic 
boundry conditions we have e{hc, eye) = L^(?>Ly — l)/2, 
where must be even. By ea. (|3.3|l . 



p{hc, eye) 



2(1, - 3) 



1 



(3.25) 



(a) 




In the thermodynamic limit, the radius of the approxi- 
mating circle decreases, approaching unity from above. 
For the shift of the circle, by eq. H3.15|l . 



s{hc, eye) = ^ 



4(3 - L,) 
3L,{3Ly - 1) 



(3.26) 



The zeros of the cyclic strip with Ly = A and m = 9{L.j; = 
18) in the x^c plane when q = 1000 are shown in Fig. 
Ol^b). The radius of this approximating circle is about 
1.42. 

Since we represent the honeycomb lattice in the form 
of bricks oriented horizontally, the results for cylindri- 
cal boundary conditions cannot be obtained from cyclic 
boundary conditions by simply switching and Ly. We 
have e{hc,eyl) — Ly{3Lx/2 — 1) for the honeycomb lat- 
tice with cylindrical boundary conditions where Ly must 
be even. When q is large, we calculate 



p{hc, eyl) = 



2(2Lj, - 3) 
3Ly[3Lx — 2) 



and 



s(/ic, eyl) 



1 4(3 - 2Ly 



3 LyiSL, - 2) 



(3.27) 



(3.28) 



The zeros of the cylindrical strip with Ly = 4 and L^ = 9 
in the Xhc plane when q — 1000 are shown in Fig. I^c). 
The radius of this approximating circle is about 1.26. 
We have also considered other boundary conditions, as 
shown in Fig. |21 



(b) 




(c) 




(d) 

FIG. 3: Zeros of the Potts model partition function in the Xf^^ 
plane for the honeycomb lattice when q = 1000, for the following 
boundary conditions and sizes; (a) toroidal with Ly = A and Lx = 
18; (b) cyclic with Ly = 4 and L^ = 18; (c) cylindrical with Ly = i 
and Lx = 9; (d) free boundary conditions with Ly = A and Lx = 9. 
For comparison, the unit circle and approximating circles are also 
shown. 



E. Kagome lattice 



As an example of a heteropolygonal two-dimensional 
lattice (an Archimedean lattice comprised of more than 
one type of regular polygon), we consider the kagome 
lattice. For a section of this lattice with toroidal bound- 
ary conditions, n{kag,tor) = SL^Ly and e{kag,tor) — 
QL^Ly so that n(kag,tor)/e{kag,tor) = 1/2 as for the 
square lattice. Hence, 

Xkag = ^ (3.29) 

and 

p{kag,tor) = - . (3.30) 

The radius of the approximating circle increases as the 
size of the lattice increases and approaches unity from be- 
low in the thermodynamic limit. The zeros of the toroidal 
strip with Ly = 2 and L^; = 9 in the Xkag plane when 
q — 1000 are shown in Fig. 01a). The radius of this 
approximating circle is about 0.938. We have also con- 
sidered other boundary conditions as shown in Fig. ^ 

IV. SUMMARY 

In summary, we have presented some new results on 
the distribution of complex-temperature zeros of the par- 
tition function of the q-state Potts model for q — > oo. 
Generalizing previous work of Wu and coworkers for 
Cartesian lattices, we have shown that with an appropri- 
ate definition of the thermodynamic limit and the limit 
q — > CX3 on an arbitrary regular lattice A of dimensional- 
ity d > 2, the zeros lie on the circle |a;A| = 1. We have 
also studied the distribution of zeros for finite sections of 
various two-dimensional lattices for large q, showing how 
they also lie approximately on circles. 
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(b) 



(d) 

FIG. 4: Zeros in the Xf^^g plane for Potts model partition function 
on the Kagome lattice when q = 1000, for the following boundary 
conditions and sizes; (a) toroidal with Ly = 2 and = 9; (b) 
g cyclic with Ly = 3 and Lx = 9; (c) cylindrical with Ly = 2 and 
Lx = 9; (d) free with Ly = 3 and Lx = 9. For comparison, the 
unit circle and approximating circles are also shown. 



V. APPENDIX I: ID LATTICE AND QUASI-ID 
LATTICE STRIPS 

Our derivation of eq. (|2.15(l assumed that the lattice 
graph G has the property that removing a few edges does 
not cut it into disconnected parts. This property does 
not hold for a ID lattice or quasi-lD lattice strips of 
sufficiently small width. Here we comment briefly on this. 
We consider first a one-dimensional lattice (line graph) 
with free or periodic boundary conditions, labelled FBC 
and PBC, respectively. We denote such graphs with n 
vertices as L^^f and L^^p- The coordination number 
K = 2 for Ln^p, and this is the effective value of the 
coordination number also for when n —^ oo, so that 
in both cases, xiu = v/q. From elementary calculations 
one has the exact results 

= + (5.1) 

and 

Z(L,,p, q, v) = <z"[(l + x,nr + ^>id] (5.2) 

As in the text, we focus on large n and large q. Aside 
from the zero at q — 0, which is not relevant in this case, 
the zeros of Z{Ln,F, q, v) occur at the single point 

xiD = -1 (FBC) . (5.3) 

Note that the number of edges for this graph is e(L„^F) = 
n — 1. Because the degree in q of the coefficient Cj{q) 
defined in eq. (|2.1(l is pj = n — j, all terms have the same 
order and none can be neglected. 

In the limit n —^ oo, followed by g — )■ oo, the zeros 
of Z{Ln^p, q, v) are given by the solution to the equation 
|1 + xid\ = |2;id|, which is the vertical line Re{xiD) — 
— 1/2 in the xid plane, i.e., 

xiD = -^+W (PBC). (5.4) 

where —oo<y<oo. As these exact results show, the 
ID case is different from the thermodynamic limit of lat- 
tices with dimensionality d > 2 in that the locus of zeros 
is strongly dependent upon the choice of boundary con- 
ditions and, moreover, is not the circle \xa\ — 1. 



In contrast, for strips of the square lattice with various 
types of boundary conditions that maintain the duality 
of the infinite square lattice, we have shown via exact 
solutions in Ref. that as g oo, the zeros do lie 
exactly on the unit circle |a;A| — 1. For these strips we 
also found that for large but finite q, the locus B consists 
of the union of (i) a self-conjugate portion of this unit 
circle given by Xgq = e*^ with Oq < 6 < tt and — tt < 6* < 
-~9o (where the angle depends on the strip) with (ii) a 
line segment —p < x^q < —1/p on the negative real axis, 
where p > 1 is a positive real constant (see the plots 
given in Ref. [l^). The arc of the unit circle and the 
line segment intersect at Xsq = —1. As q oo, p —^ 1 
and Oq ^ 0. Thus, for large but finite q, B is almost 
the unit circle except for a small line segment centered at 
Xgq = — 1 and a gap in the circle in the vicinity of Xgq = 1. 
As we discussed in Ref. these deviations can be 

understood from the fact that the large-g expansion of 
the dominant eigenvalues of the transfer matrix contain 
poles at Xsq = ±1, showing that regardless of how large q 
is, the large-g expansion breaks down at these two points. 

More generally, we have carried out similar large-g ex- 
pansions for the dominant eigenvalues of the transfer ma- 
trix for quasi-one-dimensional strips of several lattices 
with different boundary conditions. We find that for the 
strips we have considered, these expansions exhibit sin- 
gularities on the unit circle xa — 1 for strips with peri- 
odic transverse boundary conditions but not for strips 
with free transverse boundary conditions. For exam- 
ple, we find that this expansion for the Ly ^ 2 cylin- 
drical strip of the square lattice has branch-point singu- 
larities at Xsq = ±e"/^ and Xsq = ±e~"/^. Similarly, 
for the Ly — 2 cylindrical strip of the triangular lat- 
tice we find that the large-g expansion of the dominant 
eigenvalues has branch-point singularities at Xtn = ±1) 
xtri = ±e'''/^ and xtri = ie"*''/^. 

These deviations are characteristic of certain quasi- 
one-dimensional strips, which do not satisfy the premises 
of our general analytic results in Section II. We do not 
find such deviations for the thermodynamic limits of sec- 
tions of regular lattices of dimensionality d > 2 in the 
limit of large q. 
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